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The suppression of rectification at metal-Mott-insulator interfaces, which was previously 

■ shown by numerical solutions to the time-dependent Schrodinger equation and experiments 

■ on real devices, is reinvestigated theoretically using nonequilibrium Green's functions. The 
one-dimensional Hubbard model is used for a Mott insulator. The effects of attached metallic 
electrodes are incorporated into the self-energy. A scalar potential originating from work- 
function differences and satisfying the Poisson equation is added to the model. For electron 
density, we decompose it into three parts. One is obtained by integrating the local density 
of states over energy to the midpoint of the electrodes' chemical potentials. The others, 
obtained by integrating lesser Green's functions, are due to couplings with the electrodes and 

pH . correspond to an inflow and an outflow of electrons. In Mott insulators, incoming electrons 
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and holes are extended over the whole system, avoiding further accumulation of charges 
relative to that in the case without bias. This induces collective charge transport and results 
CSI ■ in the suppression of rectification. 
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^ ■ 1. Introduction 
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Correlated electron systems can be candidate materials for novel functions of electronic 
devices. Most theories for electronic devices are, however, based on a conventional one-electron 
picture for band semiconductors or band insulators. For correlated electron systems, such a 
picture may need to be replaced. Devices are always surrounded by interfaces, where the 
match or mismatch of work functions or bands^ must be considered in an appropriate man- 
ner. Generally, where two materials with different work functions are attached, bands are 
bent. A Schottky barrier is formed at a metal-insulator interface,^' ^ which is governed by the 
long-range Coulomb interaction. Because the barrier height is modulated by external bias, 
the current is usually not an odd function of external bias. The application of forward (re- 
verse) voltage lowers (raises) the interfacial barrier, leading to a larger (smaller) current, i.e., 
rectifying action at the metal-band-insulator interface.*^ 

*E-mail: kxy@ims.ac.jp 
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The importance of electron correlation in charge transport through metal-insulator in- 
terfaces is recognized in metal-insulator-semiconductor field-effect transistor (MISFET) de- 
vice structures based on organic single crystals of the quasi-one-dimensional Mott insulator 
(BEDT-TTF)(F2TCNQ) [BEDT-TTF=bis(ethylenedithio)tetrathiafulvalene, F2TCNQ=2,5- 
difluorotetracyanoquinodimethane].^ In Mott insulators, the carrier injections are ambipolar 
even if the work function of the crystal is quite different from that of the electrodes. Such char- 
acteristics are reproduced in the one-dimensional Hubbard model attached to a tight-binding 
model, where the formation of Schottky barriers is taken into account by added potentials sat- 
isfying the Poisson equation.^ A connection was suggested between these ambipolar field-effect 
characteristics and the suppression of rectification at metal-Mott-insulator interfaces.'' 

We have shown that rectification at metal-Mott-insulator interfaces is indeed sup- 
pressed, compared with rectification at metal-band-insulator interfaces, even for large work- 
function differences.^ This fact is demonstrated by numerical solutions to the time-dependent 
Schrodinger equation and also by experiments on real devices made of organic crystals of 
(BEDT-TTF)(F2TCNQ). With the current-voltage characteristics I=f{V) of the metal- 
insulator interface, it is shown that the drain current /d of the field-effect transistor is approx- 
imately proportional to VofiUc) for small Vd; Vd denotes the drain voltage, Uq the gate 
voltage in the symmetric-gate operation,^ and f'{V) is the derivative of f{V). Thus, the am- 
bipolar field-effect characteristics [V^/'(?7g) ^ even function of Uq] are the consequence 
of the suppressed rectification [/(V) ~ an odd function of V] at the metal-Mott-insulator 
interfaces. 

In the above-mentioned theoretical studies, after the ground state is obtained, a finite volt- 
age is suddenly applied and maintained at a constant value. Because we employ the periodic 
boundary condition for a finite system (with appropriate choice of a gauge) , the current under 
a constant bias finally oscillates owing to the finite-size effect.^ Thus, the current density is 
estimated by averaging the time-dependent one over a given time period. The systems that 
can be treated by this method are limited to those systems in which the band structure of the 
left electrode coincides with that of the right electrode by shifting a constant energy. Further- 
more, the time-averaged quantities are generally different from the corresponding quantities 
in the steady state. 

In order to avoid such artifacts, we need to treat steady states, which are free of the 
limitation on the metallic electrodes and the current oscillation due to the finite-seize effect. 
Without a heat bath, a steady state is reached by couplings with infinitely large, metallic elec- 
trodes. In studying the current-voltage characteristics in such a steady state, nonequilibrium 
Green's functions are common tools. Because we employ metallic electrodes where electrons 
are noninteracting, the effects of the electrodes can be incorporated into self-energies. ^'^"^^ 
They allow us to discuss the characteristics in terms of the electron density distribution. 
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In this study, we employ them to reinvestigate the current-voltage characteristics caused 
by metal-Mott-insuIator interfaces. The suppression of rectification at these interfaces can 
be elucidated by observing the spatial dependence of a "nonequilibrium" part of electron 
density. We will reformulate the Green's-function method in such a way that we distinguish 
an inflow of electrons from an outflow of electrons in the "nonequilibrium" part of electron 
density, which is crucial to the interpretation of numerical results. We will show that, in Mott 
insulators, incoming electrons and holes are extended over the whole system, giving rise to a 
collective charge transport responsible for the suppressed rectification. 

2. Model and Method 

We consider an insulator, to which metallic electrodes are attached on the left and right 
sides. For this central part, we use the one-dimensional Hubbard model (with on-site repulsion, 
U >0) for a Mott insulator and the one-dimensional tight-binding model with alternating 
transfer integrals {5t 7^0) for a band insulator, both at half filling: 

Lc 

Hcen = J2 [^^"^ + ^("^T - 1/2) (na - 1/2)] 

i=l 

Lc-1 

- E E + (-1)'<^*] + 4+l,aCi,a) , (1) 

1=1 cr 

where (cj,o-) creates (annihilates) an electron with spin a at site i, = cl^Ci^^, and 
Ui = X^cr'T'io-- Lc is the total number of sites in the central part. The scalar potential ipi is 
introduced above to account for the redistribution of electrons at interfaces to form barriers, 
compensating for the work-function differences (pi and in equilibrium.^ The applied voltage 
V is defined such that it is positive when the right electrode has a lower potential (for electrons) 
than the left and the current (without multiplication of charge) fiows to the right. 

Although the potential i/ji is defined on lattice points, we solve the Poisson equation in 
the continuum space, 

j} = -Vp{{n)-l) , (2) 

where the potential and the expectation value of the electron density per site (n) are 
functions of x, and Vp comes from the long-range Coulomb interaction. In order to match the 
Fermi levels, we set the boundary condition, i.e., the potentials in the metallic electrodes, as 

ip{x) = -(t>L + V/2 for ar < 1 , 

V'(x) = -ct>R-V/2 iovx>Lc. (3) 

In order to solve the Poisson equation analytically, we assume —dn{tp) / dip = k with a con- 
stant compressibility k, as in ref. 8 where analytic formulas are shown. There are alternative 
approaches tested, as explained in Appendix A. We have confirmed using the expectation 
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value (n) obtained self-consistently that the current-voltage characteristics and the charge 
distributions are qualitatively unchanged. 

The effects of the left and right (a=L, R) electrodes consisting of noninteracting electrons 
on the central part [eq. (1)] are described generally by retarded self-energies. In the wide-band 
limit, they are independent of energy and their matrix elements with the site indices i and j 
are given by^^ 

= -(i/2) (Tahj ^ -mjaSi,iJj,i^ , (4) 

where Sij=l for i=j and otherwise, zl=1 denotes the site connected to the left electrode, 
and iR=Lc denotes the site connected to the right electrode. Within the Hartree-Fock ap- 
proximation, the retarded Green's function for spin a, G'^{e), is given by 

[G:ie)-\ = e6.,, - , (5) 

with 

{H^ah = {HHFa)i,j ' (i/2) E ^o.Si,iJj,i^ , (6) 

a=L,R 

where the diagonal elements of Hupa are given by 

{HnFa)i,i = ^i + U{nia - 1/2) , (7) 
with a=—a, and the off-diagonal elements are 

{HHFa)i,i+i = {HnFa)i+i,i = - [U + {-lYSt] , (8) 

and {HHFa)ij=0 for | i - j \> 1. 

We numerically solve the eigenvalue equation^^ 

Lc 

5^ m^,),. u^(j) = - i7;;./2) u-^ii) , (9) 

where e^^ — 17^/2 (e^ and are real) is an eigenvalue of the complex symmetric matrix 
-^HFcr corresponding to the right eigenvector u^(i). We can set up the eigenvalue equation for 
the adjoint matrix^^ 

Lc 

J2 mFa)ij <U) = {<n + i7^/2) vl,{i) , (10) 

.?=i 

where -f^HFo- i^ Hermitian conjugate of H^p^. The u's and u's are not identical because 
^HFa -^HFo- non-Hermitian matrices. Bccaiisc H^^^ and H^^^ are symmetric matri- 
ces, we can take v'^{i)=u'^* (i). In other words, the left eigenvector is given by the complex 
conjugate of the right one. They arc functions of densities (nj^) (^=1, • • • , Lq; cr=], j) and 
must be determined self-consistently. The eigenvectors are normalized according to 

J2<*(i)<ii) =J2<ii)<(^) = ^mn . (11) 
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For the numerical solutions presented below, we can always find a complete set of eigenvectors 
satisfying 

= T.<(i')<(j') = ■ (12) 

m m 

Then, the retarded Green's function is written as 

For the density {riia-), we modify the frequently used formula in ref. 14 so as to maintain 
symmetry concerning simultaneous particle-hole transformation (cj^o- i^^Y^^la) space- 
inversion operation (i Lq + 1 — i) for ^l=^r and 4>l+4>R=0 i^C is assumed to be an even 
number). The modification is also useful for distinguishing between an inflow and an outflow 
of electrons, as will be shown later. To be more precise, we decompose the density into the 
"equilibrium" part and the parts due to the couplings with the left and right electrodes: 

a=L,R 

where the "equilibrium" part is defined by integrating the local density of states over energy to 
the midpoint of the left and right chemical potentials, fic = il^L + I^r)/'^- When the left (right) 
chemical potential is higher, (5n^ i^nf^) corresponds to the inflow, and the other corresponds 
to the outflow. The "equilibrium" part is then written as 

/^d6|-llm[G^(6)],,}/c(6) , (15) 

where /c(e) is the Fermi distribution function in a virtual system with the chemical potential 
HC- Considering zero temperature, we substitute the step function for the Fermi distribution 
function: /c(e) = ^(a*c — e)- The present decomposition is general and useful irrespective of 
whether or not the wide-band limit is applied. In general cases, the energy dependence of Fq, 
should be kept below. The definition in ref. 14 corresponds to the setting of at either 
or fXR. 

For the "nonequilibrium" part of the density, remember that the lesser self-energy is given 
in the wide-band limit by^^ 

E<{e) = i[TLfLie) + TRfR{e)] , (16) 

with /a(e) = 9{fJ-a ~ c)- It can be decomposed again into the "equilibrium" part and the parts 
due to the couplings with the left and right electrodes: 

S<(e) = S<^'i(e)+J]<5S<-(6), (17) 

a 

with 

S<^^(e) = i(ri + r«)/c(e), (18) 
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and 

5S<"(6)=ir„ [/„(€) -/c(e)] . (19) 

Because the contribution from the "equihbrium" part ^^^''(e) is regarded as included in n|^, 
the "nonequilibrium" part of the density, (5n^, can be defined as 

1 i-oo 

where the "nonequihbrium" part of the lesser Green's function, (5G^°(e), is given by the 
Keldysh equation 

5G<"(6) = G^(e)5S<"(e)G^(e), (21) 

with 6T,^"{e) in eq. (19) and G"(e) being the Hermitian conjugate of G^{e). 

Using eq. (13), we will derive each part of eq. (14). The local density of states is written 

as 

-ilmK(£)l« 

- ?'""""<^«^^(^^1TO72P- '''' 

When we substitute eq. (22) into eq. (15), the second term of eq. (22) gives a logarithmic term. 
This is an artifact of the energy-independent imaginary parts of the eigenvalues in eq. (9) (i.e., 
the wide-band limit), so that it is ignored below. The first term of eq. (22) gives 



Itan-^i^^^^ + l 



(23) 



The terms in the bracket above are reduced to the step function in the limit of 7^ ^0 (i.e.. 

After we substitute eq. (13) into eq. (20), we obtain the "nonequilibrium" part by the 
integral ranges from fj,c to /Xq- 



^ 1^ d6 [G:{e)T^GUe)U [Ue) - fc{e) 

/oo 
d6 I [Glie)]^,^ P [Ue) - fc{e)] 
-00 



7a 

27r 



7a <n{iKn(ia)K* (0<* («a) 



2^^c em - - i7^/2 - 17^/2 



X 



n,m 

1 1 



e - em + i7m/2 e - e'^ - 17-/2 

7a ] J / ''^m(^)^m('^a)'^n (0"^° 



27r I \e^-e-- 17-/2 - 17-/2 
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+ tan' 



-1 2(/ia - f-) 



7^ 



tan" 



-i2(//(--c-) 



+ Re 



X 



+ (7;^/2)' 



ilnj^ 
2 (mc 



+ (7^^./2)^ 



c 



+ (7^/2)^ 



(24) 



This expression ensures that Snf^ is a real quantity, in contrast to the formula in ref. 14. In 
the atomic hmit, eqs. (23) and (24) are simplified, as shown in Appendix B. 

Finally, by using the formula in ref. 15 with e=/i=l, the current from the left electrode is 
given by 

/oo J , 
^ ^ E [riG^(e)r«G« (6)] [/i(6) - fnie)] 



J 



ILIR 

27r 

ILIR 

27r 



J-oo „ 

/oo , 
de^{|[G^(e)]^^i |2 [/^(e)-/c(6)] 
-oo „ 



(25) 



Therefore, the current is expressed by the "nonequilibrium" parts of the density. We need to 
substitute ij,l=V/2, iiji=—V/2, and nc=0 into all the equations above. 

3. Results 

Before showing the density distributions of electrons, we compare the current-voltage 
characteristics directly obtained by time evolution through a constant potential difference 
with those obtained by the nonequilibrium Green's functions in the wide-band limit. As in 
our previous study,^ we have numerically solved the time-dependent Schrodinger equation for 
the insulator sandwiched between two metallic electrodes with different work functions and a 
common bandwidth. Figure 1(a) shows the current- voltage characteristics thus obtained for 
the Mott insulator. The parameters are introduced and the current density / is defined in ref. 8 
and Lc = L — L^. The rectifying action is shown to be suppressed. Because the bandwidths of 
the metallic electrodes are finite, the current density / tends to become saturated for a large 
V. 

Figure 1(b) shows the current- voltage characteristics obtained by the present approach 
with the nonequilibrium Green's functions for the Mott insulator. The rectifying action is 
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Fig. 1. (Color online) Current-voltage characteristics of Mott insulator with [7=1.6, 6t~0, tc^l, 
Vp=0.05, and 0^?=— 0.165 with different (pL values as indicated, (a) obtained by the time- 
dependent Hartree-Fock approximation {Lc—51, Le=49, and ie=l),^ and (b) obtained by nonequi- 
librium Green's functions (Lc=50 and 7l=7_r=1)- 



shown to be suppressed. Because the bandwidths of the metallic electrodes are assumed to 
be infinite, no saturation is observed. The values of the other parameters U, 5t, tc, Vp, 
and are the same as those in Fig. 1(a), and Lc is close to that used in Fig. 1(a). However, 
the finite metallic electrodes in Fig. 1(a) are replaced by infinite ones in Fig. 1(b), making 
direct comparison difficult. As a general trend, the current-voltage characteristics directly 
obtained by time evolution show smoother curves possibly because the current density is 
averaged over a given time period. As the system size Lc increases, the present approach 
gives smoother characteristics, as shown in Fig. 2. In any case, the suppression of rectification 
for Mott insulators is observed by both of the methods discussed above in a wide parameter 
space spanned by the coupling strength U, the system size Lc, the Coulomb parameter Vp, 
and the work- function differences <j)L and (pn. 

Figure 2 shows the current J, the "nonequilibrium" part of the density due to the coupling 
with the left electrode at the rightmost site Snj^ , and that due to the coupling with the right 
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Fig. 2. (Color online) Current-voltage characteristics, Snj^^., and (5nf , of (a) Mott insulator with 
U—1.3 and 6t—Q, and (b) band insulator with U=0 and St^O.QAd. The other parameters are 
Lc=100, tc=l, V"p=0.04, jL^jii^l, 0L=O.22, and (f)R=-Om. 



electrode at the leftmost site 5n{^, which are related by J = ^rStl^^ — ^l^Ui in eq. (25). In 
the antisymmetric case of ^l=1r and 4>l=—4'R i^ot shown), the present approach guaran- 
tees the symmetry concerning simultaneous particle-hole transformation and space-inversion 
operation, which leads to Sn'j^^_^^_-=—6n^ for any i. In a more general case of (pi 7^ —(pR but 
7l=7_R) Fig- 2 shows that the approximate relation Sn'j^^ ~ —6nf still holds for only Mott 
insulators with (pL and <j)R used here. In the band insulator, the absolute value of J is larger 
for V > [Fig. 2(b)], where the barrier is lower than that for y < 0, as expected. Later in 
Fig. 8, we will see how the rectification is realized in the band insulator. 

Prom now on, unless otherwise stated, we compare Mott and band insulators with a gap 
A=0.1, which is much smaller than the values used in previous studies. If we numerically 
solved the time-dependent Schrodinger equation as before for such a small gap, we would 
need to adopt a small V to calculate the evolution for a long time, which is proportional to 
Lq/V . Here, we use Lc=200 (larger than before) and Vp=Q.{)2> (slightly smaller than before), 
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Fig. 3. (Color online) Scalar potential ipi for T^=— 0.1 (with barriers for left-going bias) and 1^=0.1 
(without barriers for right-going bias), obtained analytically with the assumption of a constant 
compressibility. The other parameters are Lc=200, the bandwidth W^4, A=G.l, Vp=G.03, and 
<^L=-0i?=O.O5. 



which would also need long calculations in the previous approach. Because the present study 
focuses on the density distributions of incoming and outgoing electrons, we fix the parameter 
set to show them. Unless otherwise stated, the work function of the left (right) electrode is 
set to match the bottom (top) of the upper (lower) Hubbard or conduction (valence) band, 
(pL=—(j)jl=A/2. Thus, the barriers at the two interfaces almost or completely disappear for 
the right-going bias {V >0), and they become prominent for the left-going bias {V <0), as 
shown in Fig. 3. The absolute values of J are large for V >0 and small for V <0 for the band 
insulator. 

Observing the charge distribution will help to understand the mechanism for making such 
a characteristic difference between Mott and band insulators. Figure 4(a) shows the charge 
density, (nj) — l=X^o-(^«o-) — 1, and its "equilibrium" part, n^*^ — l=X^o-^iCT ~ ^ Mott 

insulator with a left-going bias and barriers. Electrons accumulate near the left interface, 
while holes accumulate near the right interface. This is because the bands are bent near 
these interfaces (Fig. 3). It should be noted that electrons and holes already accumulate for 
V=0, which is responsible for the respective band bending at the interfaces and the resultant 
matching of the chemical potentials to reach equilibrium from the isolated case with 7q.=0. By 
giving a finite bias, a certain number of electrons come in and the same number of electrons 
(holes) go out (come in) in the steady state, so that the incoming electrons or holes cause 
no instability. If one takes a closer look, the charge density alternates between even and odd 
sites, although its amplitude is very small. This 2kp oscillation is induced by the boundaries. 

Figure 4(b) shows the "nonequilibrium" parts of the charge density, one of which is due 
to the coupling with the left electrode, ^nf=Y^„^'nf^, and the other is due to that with the 
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Fig. 4. (Color online) (a) Charge density, (n^) — 1 (red or gray), and "equilibrium" part, n°'^~ 1 (light 
blue or light gray), and (b) "nonequilibrium" part due to coupling with left electrode, Snf (<0, 
light blue or light gray), and that with right electrode, Snf (>0, red or gray), for Mott insulator 
{U=l, St=0) with left-going bias (V^==-0.1). The other parameters are Lc=200, tc^l, Vp=0.03, 
7L=7fl=0.1, and 0l=-07?=O.O5. 

right electrode, Snf=Y^^ "^^f^- Here, the voltage is negative so that electrons are left-going on 
average. Because of the coupling with the left electrode, electrons go out and 6n^ <0. Owing 
to the coupling with the right electrode, electrons come in and Snf >0. Thus, the present 
definition of 5nf is useful for distinguishing between these flows. The property unique to Mott 
insulators is that 6nf is extended over the whole system in such a manner that 6nf is almost 
constant except the 2kF oscillation. Mott insulators disfavor further accumulation of charges 
relative to that in the case of V=0, so that incoming electrons and holes are extended over 
the system. 

Figure 5 shows (nj) — 1, n^^ — 1, Snf, and Snf for a Mott insulator with a right-going 
bias. For this particular set of {cpL, cpR, y)=(0.05, —0.05, 0.1), the potential ^pi is constant 
everywhere [eq. (3)] and the density (nj) is almost unity (note the vertical scale) except for 
the 2kF oscillation [eq. (2)]. Here, the voltage is positive so that electrons are right-going on 
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Fig. 5. (Color online) (a) Charge density, (n^) — 1 (red or gray), and "equilibrium" part, n°'^~ 1 (light 
blue or light gray), and (b) "nonequilibrium" part due to coupling with left electrode, Snf (>0, 
light blue or light gray), and that with right electrode, Snf (<0, red or gray), for Mott insulator 
{U=l, 5t—0) with right-going bias {V—0.1). The other parameters are the same as those in Fig. 4. 

average. Electrons come in from the left, Snf >0, and go out to the right, 6nf <0. Both 
6nf and 6nf are almost constant except for the small 2kp oscillation. In Mott insulators, 
the averaged (i.e., 2A;i?-oscillation smoothed out) densities of the incoming electrons and holes 
show not only a small spatial modulation but also insensitivity to the sign of V. This fact 
ensures the suppression of rectification because the current is given by the difference between 
the density of the incoming electrons at the exit and that of the outgoing electrons at the 
entrance [multiplied by the respective 7q,, eq. (25)]. 

The behavior of the metal-Mott-insulator interfaces is contrasted with that of the metal- 
band-insulator interfaces. Figure 6 shows (jij) — 1, n^^ — 1, Snf, and Snf for a band insulator 
with a left-going bias and barriers. It is clearly shown in n^^ — 1 that electrons accumulate 
near the left interface, while holes accumulate near the right interface [Fig. 6(a)]. They are 
much larger than the corresponding quantities at the metal-Mott-insulator interfaces. Thus, 
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Fig. 6. (Color online) (a) Charge density, (n^) — 1 (red or gray), and "equilibrium" part, n°'^~ 1 (light 
blue or light gray), and (b) "nonequilibrium" part due to coupling with left electrode, Snf (<0, 
light blue or light gray), and that with right electrode, 5n^ (>0, red or gray), for band insulator 
{U=0, (5^=0.025) with left-going bias {V=—0.1). The other parameters are the same as those in 
Fig. 4. 



the accumulation is more sensitive to the band bending. 

The difference from the metal-Mott-insulator interfaces is conspicuous in Snf and Sn^ 
[Fig. 6(b)]. Here, electrons are left-going on average. The quantity 6nf is thus negative, but 
its magnitude becomes almost zero at some point. On its left side, —Snf is largest at the 
left interface and decays with increasing distance from the interface. This behavior originates 
from the outgoing electrons. On its right side, —6nf shows a maximum away from the right 
interface. This behavior can be regarded as due to ipi, which is higher on the right side. In other 
words, holes further accumulate on the right side at a negative V. The spatial dependence of 
dnf is obtained by the particle- hole transformation and the space- inversion operation of Snf. 
The quantity Snf is positive and becomes almost zero at another point. On its right side, Snf" 
is largest at the right interface and decays with increasing distance from the interface. This 
behavior originates from the incoming electrons. On its left side, Snf shows a maximum away 
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Fig. 7. (Color online) (a) Charge density, (n^) — 1 (red or gray), and "equilibrium" part, n°'^~ 1 (light 
blue or light gray), and (b) "nonequilibrium" part due to coupling with left electrode, Snf (>0, 
light blue or light gray), and that with right electrode, 5nf (<0, red or gray), for band insulator 
{U=Q, (5i=0.025) with right-going bias (V^=0.1). The other parameters are the same as those in 
Fig. 4. 



from the left interface owing to further accumulation of electrons at lower 'i/'i's- This further 
accumulation is realized by the absence of on-site repulsion. 

Figure 7 shows (rij) — 1, nj^^ — 1, 6nf, and 6n^ for a band insulator with a right-going bias. 
For this particular set of {(pL, 4>Ri ^)) the potential -01 is constant. The quantity rf^ is unity, 
but (rej) largely deviates from it. Because the density of electrons is much more sensitive to 
the bias-induced change in the potential distribution for band insulators, (nj) and n^^ are 
now quite different. Here, electrons are right-going on average, so that 5nf >0 and Snf <0. 
Because of the absence of barriers, a substantial number of electrons come in and go out. 

In order to see how the rectification is realized in a band insulator, we show its 5nf and 
5nf in Fig. 8 with parameters used in Fig. 2(b) and different biases, whose magnitudes are 
larger than the gap A=0.18 here. For V < —A, the absolute value of J, | J |, is smaller than 
that for V > A [Fig. 2(b)]. In this case, the density of the incoming electrons 6nf > decays 
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Fig. 8. (Color online) "Nonequilibrium" part of density for (a) < {Snf- for y=— 0.24, Snf^ for 
y=-0.20, Snf for F=-0.20, and 6nf for y=-0.24 from top to bottom) and (b) > {Snf 
for V=0.2A, Snf for y=0.20, Snf for y=0.20, and Snf for y=0.24 from top to bottom) in band 
insulator with U=0 and 5t~Q.045. The other parameters are the same as those in Fig. 2. 



with increasing distance from the right interface, while the density of the outgoing electrons 
Snf < decays with increasing distance from the left interface. With increasing | ^ |, both 
I Snf I and | Snf \ slightly increase at each site i, but they decay in quite similar manners 
[Fig. 8(a)]. 

For y > A, on the other hand, | J \ increases more rapidly with | V \ than for V < —A. 
Now, the current substantially flows through the system. The density of the outgoing electrons 
Snf < decays with increasing distance from the right interface. However, with increasing 
I y I , a substantial number of electrons are shown to penetrate from the left interface into the 
system [Fig. 8(b)]. Both | Snf \ and | Snf \ finally decay toward the right and left interfaces, 
respectively, giving small | Snf^^ \ and | Snf \ values compared with | Snf \ and | Snf^ |, 
respectively. However, | Snf^ \ for y > A is much larger than | Snf^ \ and | Snf \ for 
V < —A because of the penetration of a larger number of electrons into the system. These 
behaviors of Snf and Snf in a band insulator are in contrast to those in a Mott insulator. 
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where the averaged (i.e., 2A;ir-oscination smoothed out) densities of the incoming electrons 
and holes show not only little spatial modulation but also insensitivity to the sign of V. The 
V dependence of 6nf and that of 6nf in the Mott insulator closely follow that of 5n£^ and 
that of 6nf- in Fig. 2(a), respectively. 

4. Summary and Conclusions 

To elucidate the suppression of rectification at metal-Mott-insulator interfaces, which 
is obtained by numerical solutions to the time-dependent Schrodinger equation in our pre- 
vious study,^ we employ nonequilibrium Green's functions in this paper. We consider one- 
dimensional half-filled electron systems and use the mean-field Hubbard model for a Mott 
insulator and the transfer-alternating tight-binding model for a band insulator. Metallic elec- 
trodes consisting of noninteracting electrons are attached to the system, and their effects are 
incorporated into self-energies within the wide-band limit. We take account of work-function 
differences, which are responsible for the formation of Schottky barriers, by adding to the 
model a scalar potential that satisfies the Poisson equation. 

In the mean-field approximation, the retarded and advanced Green's functions are ob- 
tained by solving the eigenvalue equation for a complex symmetric matrix, whose imaginary 
part comes from the self-energies due to the couplings with the electrodes. For the electron 
density to be determined self-consistently, we modify the frequently used formula in ref. 14 so 
as to maintain some symmetry in the particular case used here. The modification is useful for 
distinguishing between an inflow and an outflow in the "nonequilibrium" part of the density. 
The current is given by the difference between these flows, which are measured at appropriate 
sites and multiplied by couplings with electrodes. 

By plotting the spatial dependence of the density of incoming electrons and holes, we 
clarify the difference between metal-Mott-insulator and metal-band-insulator interfaces. In 
Mott insulators, the incoming electrons and holes are extended over the whole system in 
such a manner that Snf and Snf are almost constant over i=l, • • • , Lc- Electrons/holes do 
not accumulate at any place, so that the bias-induced change in the density of electrons is 
extended. Thus, charge transport becomes collective. 
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Appendix A: Alternative Approach to Band Bending 

As is well known in electromagnetics, the Poisson equation is equivalent to the long-range 
Coulomb interaction, ^j^i{V/ \ rj — |) {rij — 1) (n^ — 1), where is the position vector for 
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site i. Then, we can define the Hartree potential as 




((%■) - 1) , 



(A-l) 



i^i = ipi + ai + h , 



(A-2) 



where the constants a and h are determined so that the potential V' satisfies the boundary 
condition [eq. (3)]. In this way, we can treat the band bending at interfaces and the effect of 
the long-range Coulomb interaction on the insulator (e.g., charge ordering if it is strong) on 
the same footing. This method can be extended to include the screening effect. In addition, 
we found that its numerical convergence is generally better than that of solving the Poisson 
equation. 

Appendix B: Densities in the Atomic Limit 

In the limit oitc + {—iy5t ^0, i.e., in the atomic limit, one can simplify the expressions for 
the density. The matrix i^HFo- is diagonal and its elements are given by e?^ = tpi + U [ni^ — 1/2), 
so that the eigenvalues — i7^/2 are written as e^^ — i'^a/'^ for m=ia and as otherwise. 
Then, the "equilibrium" and "nonequilibrium" parts of the density at i=ia are then written 
as 




(B-1) 



and 



6nf _ = — tan" 

a ^ 

The total density is thus written as 



i2(/^a-eU 



— tan 



-1 



2(Mc-efJl 



(B-2) 



which is a reasonable expression in the atomic limit. 




(B-3) 
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